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Abstract 



, Two new concepts, generic regular decomposition and regular-decomposition-unstable 

CN ■ (RDU) variety for generic zero-dimensional systems, are introduced in this paper and an al- 

^) ' gorithm is proposed for computing a generic regular decomposition and the associated RDU 

D , variety of a given generic zero-dimensional system simultaneously. The solutions of the given 

^0 1 system can be expressed by finitely many zero-dimensional regular chains if the parameter 

value is not on the RDU variety. The so called weakly relatively simplicial decomposition 
plays a crucial role in the algorithm, which is based on the theories of subresultant chains. 
Furthermore, the algorithm can be naturally adopted to compute a non-redundant Wu's 
Q^) | decomposition and the decomposition is stable at any parameter value that is not on the 

■ RDU variety. The algorithm has been implemented with Maple 15 and experimented with 

^ 1 a number of benchmarks from the literature. Empirical results are also presented to show 

, the good performance of the algorithm. 



> 

(N 



O 

(N 



Keywords: generic zero-dimensional system, regular-decomposition-unstable variety, para- 
metric triangular decomposition, generic regular decomposition 



1 Introduction 

Solving parametric polynomial systems is usually a key problem in many research and applied 
areas, such as automated geometry theorem deduction, stability analysis of dynamical systems, 
robotics and so on J2JE2HI]. By "solving", we often mean to determine (1) for what parameter 
values the polynomial system has solutions, and (2) whether the solutions can be expressed by 
some simple representations. 

1 Generally speaking, there are two kinds of methods for solving the above questions (1) and 

(2), i.e., the methods based on Grobner bases [101 1131 1141 121] and triangular decompositions 

mniaianniniiisEziiHiisEiEfl- 

For parametric systems, the concepts of comprehensive Grobner system (CGS) and com- 
prehensive Grobner bases (CGB) introduced by Wcispfcnning in [2T] and the algorithms for 
computing them (TDl [13j [14l [15l [THJ [2TJ are powerful tools for answering questions (1) and (2). 
The first CGB algorithm introduced in [2T] suffers from the problem of too many redundant 
branches. Many improved algorithms have been proposed since then [TU] EH HH [13 US] , among 
which, the one proposed by Suzuki and Sato [T5] was accepted widely by subsequent researchers. 
The latest progress on this subject was reported by Kapur et al. [TU]. They solved the fa- 
mous P3P problem [7] by computing CGS and provided empirical data illustrating that the CGS 
method could solve practical problems in amazingly short time. 

The methods based on triangular decompositions have been studied by many researchers 
since Wu's work [33] . A significant concept in the theories of triangular sets is "regular chain" 
(or "normal chain" ) introduced by Kalkbrener [9] and Yang and Zhang [27] independently. Gao 
and Chou proposed a method in [6] for identifying all parametric values for which a given system 
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has solutions and giving the solutions by p— chain^j] without a partition of the parameter space. 
Wang generalized the concept of regular chain to regular system and gave an efficient algorithm 
for computing it [TBI EH HHj ■ It should be noticed that, due to their strong projection property, 
the regular systems or series computed by RegSei@ may also be used as representations for 
parametric systems. The concept of comprehensive triangular decomposition (CTD) introduced 
by Chen et al. in [3] can answer questions (1) and (2). Algorithms for computing regular chain 
decompositions and CTDs have been implemented as central functions of RegularChains library 
in Maple 15. 

For a given parametric system P with n variables and d parameters, many existing algorithms 
for computing regular decomposition over a certain ring K will give a regular zero-decomposition 
of P in K Then, if one wants to answer questions (1) and (2), one may try computing 
projections from the solution space (the zero-decomposition) to the parametric space by the 
methods introduced in [U[T5]. On the other hand, there are some other methods, such as Wu's 
method [22] and relatively simplicial decomposition (RSD) [27], which consider parameters as 
"constants" during the process of decomposition and can obtain zero-decompositions of P in 
K[U] where U stands for the d parameters. In this paper, we follow the idea of the latter 
methods and propose an algorithm for computing a so-called generic regular decomposition T 
of a generic zero-dimensional system P in K[U] (see Definition 0]) . At the same time, the 
algorithm will obtain a parametric polynomial such that the regular decomposition is stable 
at any parametric point outside the variety generated by the parametric polynomial and we 
call the variety regular-decomposition-unstable (RDU) variety. Roughly speaking, "stable at a 
parametric point" means that the regular decomposition will keep after we substitute the point 
for the parameters in P and T (see Definition 0]) . As a result, questions (1) and (2) for generic 
zero-dimensional systems are answered except for the case when parameters are on the RDU 
variety. That's why the decomposition is called generic regular chain decomposition. 

The proposed algorithm is based on weakly relatively simplicial decomposition, a new concept 
that is weaker than relatively simplicial decomposition proposed by Yang et al. in |27| and in- 
spired by the method for computing regular systems introduced by Wang in [18LI19] . In addition, 
the proposed algorithm can be naturally adopted to compute a non-redundant Wu's decompo- 
sition for a given generic zero-dimensional system. Furthermore, computing RDU varieties can 
be regarded as the first step of computing border polynomial (BP), which is a crucial concept 
introduced by Yang et al. [5U (55[ [55] for solving the real root classification (RRC) problem 
of parametric semi- algebraic systems. As a matter of fact, a RDU variety of a generic zero- 
dimensional system with respect to (w.r.t.) a generic regular decomposition is a subvariety of the 
hypcrsurface generated by a certain BP. The new algorithm has been implemented on the basis 
of DISCOVERER [53J with Maple 15 and experimented with a number of benchmarks from the 
literature [U [TUl [Tj2 Q3] . Empirical results are also presented to show the good performance 
of the algorithm. 

The paper is organized as follows. Section [2] gives basic definitions and concepts that are 
needed to understand the main algorithm. Section [3] contains the main algorithm, namely Al- 
gorithm [3J and some relative subalgorithms, especially the subalgorithm for computing weakly 
relatively simplicial decompositions. Besides, proofs for these algorithms are presented in this 
section and an illustrative example is given. The empirical data and comparison with previous 
works along with several implementation details are presented in Section [4] Section 5 concludes 
the paper with a discussion on our future work along this direction. 

2 Preliminaries 

The following paragraphs give a brief outline of the vocabulary and tools we will be using 
throughout the paper. All concepts without precise definitions can be found in J2[[22l[26]. Several 
illustrative examples are presented. In these specifical examples, R and C stand for the feild of 
real numbers and the field of complex numbers, respectively. 

Suppose {iti, . . . , Ud, Xi, . . . , x n } is a set of indeterminates with a given order u\ -< ... -< Ud -< 
x\ -< ... -< x n where {ui, . . . ,Ud} and {x\, . . . , x n } are the sets of parameters and variables, 

1 The concept p— chain is stronger than regular chain, see more details in [6]. 
2 http: / /www-calfor. lip6.fr/~wang/epsilon/ 
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respectively. Let U = {u\, . . . , Ud} and X = {xi, . . . , x n }. Suppose A is a unique factorization 
domain and K is its algebraic closure. Let K[U] be the ring of polynomials in U with coefficients 
in K and let K[U] be the algebraic closure of K[U]. Let A[[/][A] be the ring of polynomials 
in X with coefficients in the ring K[U]. Let K[X] be the ring of polynomials in X with co- 
efficients in the field A. A non-empty finite subset P of A[£/][A] is said to be a system. If 
P C A[LT][A']\A[A'], it is a parametric system. If P C A[A], it is a constant system. 

For a non-empty finite subset P CA [U][X] (A [X]), (P)k[u}[x] ((P)~k [x]) denotes the ideal 
generated by P in A[[/][A] (A[A]) and sqrt((P) k[u][x]) ( sc l r t((P)7f [x])) denotes the radical ideal 
of (P) K[U][X ] «P)k[x])- For an y F in K [ U ][ X ]\{°} (K[X}\{0}) and for any x G X, if x appears 
in F, F can be regarded as a univariate polynomial in x, namely F = Cqx" 1 + Cix" 1 ^ 1 + . . . + C m 
where Co, Ci, . . . , C rn are polynomials in A [U] [X\{x}\ {K[X\{x}}) and Co 7^ 0. Then to is the 
leading degree of A w.r.t. a; and is denoted by deg(A, x). Note that if x doesn't appear in F, 
deg(A, x) = 0. If there exists p (1 < p < n) such that deg(A, x p ) > and for every i (p < i < n), 
deg(A, x^ — 0, then the class of F is p. If deg(A, a;,-) = for every i (1 < i < n), then the class of 
F is 0. The class of F in A [U] [X}\{0} (K[X]\{0}) is denoted by cls F . If cls F > 0, x chp is the 
main variable of F and is denoted by mvar(A). Assume that F = Cox™ + Cix™^ 1 + . . . + C m 
where p = c\sf > and Co 7^ 0, then Co, denoted by \p, is the initial of F and XT' , denoted by 
rank(F), is the rank of F. 

A non-empty finite set T = {T u T 2 . . . , T r } of polynomials in A"[C/][A] (K[X]) is a tri- 
angular set in A[J7][A] (A [A]) if < cls Fl < c1st- 2 < ... < clsT r . For a triangular set 
T in K[U][X] (A [A]), I T , mvar(T) and rank(T) denote n TeT lT, r {mvar(T)|T e T} and 
{rank(T)|T e T}, respectively. The saturated ideal of a triangular set T in A'[J7][A] is de- 
fined as the set {F e A[[/][A]|It s A 6 (T)k[u][x] f° r some positive integer s} and is denoted 
by sat(T) it -r £/ ir x -i. Similarly, the saturated ideal of a triangular set T in A[A] is defined as the 
set {F G A[A]|I T S F G (T)^ [Jf] for some positive integer s} and is denoted by sat(T)^r x , . 
Suppose F G A[C/][A] (A [A]) and T is a triangular set in K[U][X] (A [A]), then F is re- 
duced w.r.t. T if dcg(A, mvar(T!;)) < deg(Ti, mvar(Ti)) for every i (1 < i < r). A triangular 
set T = {Ti,T2 . . . ,T r } in A[J7][A] (A[A]) is a non- contradictory ascending chain in A[t/][A] 
(A [A]) if Ti is reduced w.r.t {Ti, . . . , Tj_i} for every i (2 < i < r). A single-element set 
{F} C A[[/] ({A} C A) is a contradictory ascending chain in K[U][X] (A [A]) if A ^ 0. Re- 
mark that an ascending chain is either a non-contradictory ascending chain or a contradictory 
ascending chain. 

For two polynomials A and P in A[[/][A] (A[A]) and a variable x £ A, the pseudo re- 
mainder |27j of A pseudo- divided by P w.r.t. x is denoted by prem(A, P, x). Particularly, 
prem(A, P, mvar(P)) is denoted by prem(A, P). For a polynomial A e A"[t/][A] (A[A]) and 
a triangular set T = {Ti, ...,T r } in A[[/][A] (A[A]), the successive pseudo remainder |27| of A 
w.r.t. T is denoted by prem(A, T), namely 

prem(A, T) = prcm(. . . prcm(prem(A, T r ), T r _i), . . . , T±). 

For a finite set P C K[U][X] (A[A]), prem(P,T) denotes the set {prem(A, T) | A G P}. 

Now we'd like to present some results introduced in [22]. For P C A[J7][A], suppose V^^(P) 
denotes the set 

{(a u ...,a„) G Kpf l \P(U,a u ...,a n ) =0,VPgP}. 

An ascending chain C in A[J7][A] is a characteristic set of P in A[J7][A] if C C (P)k\u][x] 
and prem(P, C) = {0}. According to Wu's method introduced in [221 ' f° r an y non-empty finite 
subset P C A[J7][A], there exists a finite sequence of ascending chains Ci, C2, ■ ■ • , C m (m > 1) 
in A'[t/][A] such that 

(1) C m is a contradictory ascending chain and, if to > 1, Ci, C2, ■ ■ ■ , C m _i is a finite sequence 
of non-contradictory ascending chains; 

(2) Ci, C2, ■ ■ • , C m is a finite sequence of characteristic sets in A[[/][A]; 

(3) If to = 1, V m (P) = 0. Otherwise, V^P) = ^V^C^Ic,) . 

The set of ascending chains {Ci, C2, • ■ • , C m } above is said to be a Wu's decomposition or 
characteristic set decomposition of P in A[[/][A]. In addition, if m = I or for every i (1 < 
i < to — I), mvar(Ci) = A, P is a generic zero- dimensional system. Remark that a Wu's 
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decomposition may suffer from the redundant branches problem. That means, V^ryj(Ci\Ic i ) 
can be an empty set for some non-contradictory ascending chain Gj (1 < i < m — 1) if m > 1. 

Another important concept in the theories of triangular decompositions is regular chain. For 
two polynomials F and P in A[P][A] (A'[A]) and a variable x G X, the resultant [27] of F and 
P w.r.t. x is denoted by res(P, P, 2). Particularly, res(P, P, mvar(P)) is denoted by res(P, P). 
For a polynomial F G K[U][X] (K[X]) and a triangular set T = {T u ...,T r } in K[U][X] (K[X]), 
the successive resultant [37] of F w.r.t. T is denoted by res(P, T), namely 

res(P,T) = res(. . . res(res(P, T r ), T r _i), . . . ,T 1 ). 

A triangular set T = {Ti, . . . ,T r } in A'[P][A] (A[X]) is said to be a n^uZar c/iain in A[C/][X] 
(K[XJ), if I Tl _^ and for each i (1< i < r), res(I Tj , . . . , Ti}) ^ 0. If T is a regular chain 

in A"[P][A] (A[A]) and mvar(T) = X, T is a zero- dimensional regular chain. Regular chains 
have a series of good properties, some of which are listed below. For P C AT[A], suppose V(P) 
denotes the set {(ai, . . . , a„) G K \P(ai, . . . , a„) = 0, VP G P}. 

Proposition 1. □ E H IH E! H3 HI if t is a regular chain in K[U][X] (K[X]), then Vj^T\ 
I T )^0 (V(T\I t )^0). 

Proposition 2. El El GD El EU Ell EU If T is a regular chain in K[U][X] and P is a polynomial in 
K[U][X], then 

(1) prem(P, T) = if and only if P G satfT)^^^; 

( 2 ) v kP]( T \ It ) C V l^]( p ) l / anrf ^ P e s q rt ( sat ( T )if[r/][x])- 
Furthermore, if T is zero-dimensional, then 

( 3 ) V 7f[Z7]( T ) n V iqc7]( p ) ^ */ and onl V T) = 0. 

Remark 1. El El ED EH EU EH EH This remark is an analogue of Proposition^ If ' T is a regular 
chain in K[X] and P is a polynomial in K[X], then 

(1) prem(P, T) = if and only if P G sa^TW^; 

(2) V(T\I T ) C V(P) if and only if P G sqrt(sat(T) s = [x] ). 
Furthermore, if T is zero-dimensional, then 

(3) V(T) n V(P) ^j/ and onfy i/ res(P, T) = 0. 

Remark 2. There exist various efficient algorithms for computing regular chain decompositions^ 
[7] fPl [TSl \19. 27. 28l. Regular chain decompositions don't suffer from redundant problem as Wu's 
decompositions owing to the fact that a regular chain always has solutions as shown in Proposi- 
tion [I] It should be noted that the definition of triangular set and thus that of regular chain in 
K[U][X] introduced above is not exactly the same as that introduced in J3] ^ \18\j when dealing 
with parametric systems. For example, consider a parameter system {^,21,22} in K[u] [21, 22] ■ 
The system itself is a regular chain in M.[u] [x\, X2] = R[«, xi, 22] according to the definition of 
regular chain introduced in £2 [^J \18jl and it is seen that the system has solutions if and only 
if the parameter value of u equals to zero. But {u, 21,22} is not a regular chain in WL[u][xi, X2] 
in this paper. Because we demand that mvar(T) C X for any triangular set T in K[U][X] as 
defined before and parameters Ux,...,Ud are regarded as "constants" here. This idea is inspired 
from the point of view proposed in ]22\j and plays a key role throughout this paper. We need to 
give some more precise definitions at first, then we'll turn back to show how to solve this simple 
example by our method proposed in this paper. 

Definition 1. Suppose P is a generic zero- dimensional system in K\U][X]. A finite set T of 
triangular sets in K\U][X] is said to be a parametric triangular decomposition o/P in K\U}[X], 
*/ V i?p]( p ) = Ut6tV^j<T\It). I/T = or V^(T\I T ) ^ for any T G T. the parametric 
triangular decomposition is said to be non-redundant. If T is a finite set of regular chains in 
K[U][X], the parametric triangular decomposition is said to be a parametric regular decomposi- 
tion. 

For each a = (01, . . . ,ad) G K d , <j> a : K[U][X] — > K[X] is a homomorphism such that 
4> a (F) = F(a,X) for all F G A[P][A] and we denote 4> a {F) by F(a). For a non-empty finite set 
P C K[U][X], P(a) denotes the set {F(a)\F G P} and remark that P(a) = if P = 0. 
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Definition 2. Let T be a parametric triangular decomposition of a given generic zero- dimensional 

system P in K[U][X]. For any a G it* , if V(P(a)) = U TeT V(T(a)\I T (a)) and rank(T) = 
rank(T(a)) for any TgT, then T is said to be stable at a. 

Definition 3. Let T foe a regular chain in K[U][X] and a G if . //T(a) is a regular chain 
in K [X] and rank(T(a))= rank(T), then we say that the regular chain T specializes well at a. 

Suppose V is an affine variety in K d . Then dim(V) denotes the dimension of V. Please see 
the precise definition of dimension of affine variety in [2] . 

Definition 4. Let T be a parametric regular decomposition of a given generic zero- dimensional 
system P in K[U][X]. Suppose V is an affine variety in K with dim(V) < d. If for any a G 

K d \V, V(P(a)) = U T eTV(T(o)\I T (o)) and T specializes well at a for any T G T if T ^ 0, then 
T is said to be a generic regular decomposition of P and V is said to be a regular-decomposition- 
unstable (RDU) variety o/P w.r.t. T. 

Now consider the example {u, x\, X2} in Remark[5J it is clear that is a parametric triangular 
decomposition and also a parametric regular decomposition of {u, X\,X2] in R[it] [x\, X2]. This 
parametric regular decomposition is stable at any a G C if a ^ 0. Thus is a generic regular 
decomposition of {u, x\, X2} and the variety generated by u in C is a RDU variety of {it, Xi, X2} 
w.r.t. 0. That means for almost all parameter values except some special ones, the given system 
{u,xi,X2} has no solutions. In most cases, we don't care about those special parameter values. 
But if we want to know what would happen if the parameter value equals to zero in this example, 
we can regard u, x\ and X2 as x' ll x' 2 and x' 3 , respectively. Then {2^, x' 2 , x' 3 } is a constant system 
in M[x^, x' 2 , x' 3 ] and a constant system is easy to solve. This example shows our main idea in this 
paper. In the next section, we propose an algorithm to compute generic regular decompositions 
and the associated RDU varieties simultaneously for generic zero-dimensional systems. 

To present the main result in the next section, we need to introduce some more notations. For 
any P C K{U]{X], Vj^(P) denotes the set {(a 1: . . . ,a d+n ) G 'K d+n \P(a 1 , . . . , a d+n ) = 0, VP G P}. 

For any B C K[U], V u '(B) denotes the set {(oi, . . . ,a d ) G K d \B(a u . . . , a d ) = 0, V£? G B}. For 
any F G K[U][X], the coefficients B u . . . , B t of F in X are polynomials in K[U], then V u (F) 
denotes V u ({B ll B t }). Note that for two finite subsets P and H of K[U][X], V 7? ^(P\H) 

denotes the set V^-(P)\V^-(H). Similarly, we can have V(P\H), Vk(P\H) and V l/ (P\H). 
The following Lemma[T]is proposed in Remark that the definition of regular chain in /£"[{/] [JT] 
in this paper is not exactly the same as that in [4] as mentioned in Remark[2j Therefore, Lemma 
[T] here is stated in our way. 

Lemma 1. ^ Let T be a regular chain in K[U][X]. Then T specializes well at a if and only if 
aGirV u (res(I T ,T)). 

3 Theory and Algorithm 

3.1 Weakly Relatively Simplicial Decomposition 

In this section, we introduce weakly relatively simplicial decomposition in zero-dimensional 
case, which is a weaker concept comparing to relatively simplicial decomposition proposed in 

E3- 

Definition 5. Let T be a zero-dimensional regular chain in K[U] [X] and P G K[U] [X] . Suppose 
H and G are two finite sets of zero- dimensional regular chains in K[U] [X] . If V k[u] (T U {P}) = 
UhghV^pj^H) and V^jj(T\P) — UGeGVjjpjy(G), then (H, G) is said to be a weakly relatively 
simplicial decomposition of T w.r.t. P in K\U][X]. 

Definition 6. Suppose (H, G) is a weakly relatively simplicial decomposition of a zero- dimensional 
regular T w.r.t. a polynomial P in K\U][X\. For any a G K d , if 
(1) T specializes well at a, 



5 



(2) V(T(a) U {P(a)» = U H£ iV(H(a)) and ifB. ^ 0, H specializes well at a for any Hgl, 

(3) V(T(a)\P(a)) = UGecV(G(a)) and ifG^%, G specializes well at a for any G £ G, 
then the weakly relatively simplicial decomposition (M, G) o/T w.r.t. P is said to be stable at a. 

Remark 3. A stronger concept, relatively simplicial decomposition, was firstly introduced by 
Yang and Zhang in \21\ 281. The precise definition and the algorithm can be seen in \2b\ I^Pf . 
Note that a relatively simplicial decomposition is a weakly relatively simplicial decomposition but 
the converse is not true. For instance, ({{x 2 , y}} , {{x + 1,2/}}) is a weakly relatively simplicial 
decomposition of {(x + l)x 2 , y} w.r.t. x + y but it is not a relatively simplicial decomposition of 
{(x + l)x 2 , y} w.r.t. x + y because prem(x + y, {x 2 , y}) = x ^ 0. 

Now we present Algorithm [T] for computing weakly relatively simplicial decomposition^], 
which is different from Algorithm RSD proposed in |27j . More specifically, for any zero-dimensional 
regular chain T in A'[J7][A] and for any polynomial P in A[[/][A], it is noted that either the 
first output H or the second output G of WRSD(T, P, X) can be an empty set, but MUG can't 
be empty since V^jjjj-(T) ^ 0. According to the specification of Algorithm Q] and Definition [5] 
(H, G) is a weakly relatively simplicial decomposition of T w.r.t. P. It is also noted that the 
third output F of WRSD(T, P, X) is a polynomial in K[U] and thus V u (F) is a hypersurface in 
K . The specification of Algorithm [1] indicates that for any a £ K \V U (F) , (H, G) is stable at 
a. 

In addition, we need to interpret some notations used in our pseudocodes. Assume that Alg 
is a name of an algorithm and p\ , . . . , pt is a sequence of inputs of this algorithm. If the output of 
Alg(pi, . . . ,pt) is a finite sequence q±, . . . , q s , qi is denoted by Alg(pi, . . . ,pt)i for any i (1 < i < s) 
and also said to be the i-ih output of Alg(pi, . . . ,p t ). Given a finite set S = {si, . . . , St} and a 
map on S, op(5) denotes the finite sequence Si, . . . , s t and map(s — > cf)(s),S) denotes the set 
cf>(S). 

Before showing the termination and the correctness of Algorithm [TJ we need to prepare some 
statements. In the following discussion, we assume that the readers arc familiar with the theories 
of subrcsultants. The precise definitions of subrcsultant chain and regular subresultant chain can 
be seen in [T^J IE] and Lemma [5] can be found in [5J [55] ■ 

Lemma 2. EES! Let (j) : R -> R be a ring homomorphism. Denote also by (f> the induced 
homomorphism <p : R[x] — > R[x], where both R and R are integral domains. Suppose F and G 
are polynomials in K[x] and b and c are the leading coefficients of F and G respectively. Assume 
that m = deg(F, x) > I = dcg(G, x) > and rh = deg(0(F), x) > I = deg(^>(G), x) > 0. If m > I, 
let fx = m— 1, otherwise, Jl — rh. Suppose Sj is the j-th subresultant of F and G w.r.t. x and Sj 
is the j-th subresultant of 4>{F) and 4>{G) w.r.t. x. Then 4>{Sj) — 5 ■ Sj for any j (0 < j < Jl), 
where 




4>{b) ^ and 0(c) = 0, 
W-fi^c)™-™, (f>(b) = and 0(c) ^ 0, 
<f>(b) = 0(c) = 0. 

Furthermore, if m > I, let fj, = m — 1, otherwise, let /i = m. Suppose Rj is the j-th principal 
subresultant coefficient of F and G w.r.t. x. Then (f>(G) = if 0(6) ^ and 4>{Rj) = for any 
3 (0<J<M)- 

Roughly speaking, Algorithm [1] is based on Lemma G2 which is inspired by the analogous 
results presented in JTEJ HH] • We give a proof here because the results showed in Lemma [3] is not 
covered by what is proposed in [18l [19] . Anyway, the idea of the proof is also borrowed from 

nana- 

Lemma 3. Given two polynomials F and G in K[U] [X] (0 < deg(G, x„) < deg(F, x n )), suppose 
Sd v , • ■ ■ , Sdn Sd is the regular subresultant chain of F and G w.r.t. x n . Let Sd v+1 = F. Assume 
that Rdi is the di-th principal subresultant coefficient of F and G w.r.t. x n for any i (0 < i < v+1) 
and Qdi is the pseudoquotient of F and Sd t w.r.t. x n for any i (1 < i < v). Then 



3 Lines 2 and 3 of Algorithm ^ can be removed without loss of correctness. 
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Algorithm 1. WRSD 



Input: A zero-dimensional regular chain T = {T\, . . . ,T n } in if[C/][X], a polynomial 

P £ K[U][X], variables X = {x 1; . . .,x n } 
Output: Two finite sets of zero-dimensional regular chains H and G in K [U][X] and a 

polynomial F 6 K[U], such that 
(ljVjcfoKT U {P}) = UH e HV STUT (H); (2)V^ mT (T\P) = Ug^V^G); 



(3)for any a G K \V U (F), T specializes well at a and 

(3a)V(T(a) U P(a)) = U H eH[V(H(a)) and H specializes well at a for any H e H if : 
(3b)V(T(a)\P(a)) = U GeG V(G(a)) and G specializes well at a for any G e G if G 

1 H:=0, G:=0, F:=rcs(I T ,T) 

2 if P is not reduced w.r.t. T then 

3 |_ return WRSD(T, prem(P, T),X) 

4 if P = then 

5 |_ return {T}, 0, F 

6 if clsp = then 

7 |_ return 0, {T}, P ■ F 

8 if clsp 7^ n then 
W:=WRSD({Ti, • . • , T dsp }, P, { Xl , . . . , x clsp }) 
H:=map(t i U {T clsp+1 , . . . , T„}, Wi) 
G:=map(i ->■ i U {T clsp+1 , . . . , T„}, VF 2 ) 
return H, G, F • VK 3 

13 if res(P, T) ^ then 

14 | F:=F • res(P, T), G:={T} 
is else 



^0 



10 

11 

12 



16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 

35 
36 
37 
38 
39 
40 



compute the regular subresultant chain Sd„, ■ • ■ , Sd 1 , Sd a of T n and P w.r.t. x n 
if n = 1 then 
H:={{S dl }} 

Let Q be the pscudoquoticnt of T\ and S^i w.r.t. x\ 
G:=WRSD({g}, P, X) 2 , F:=WRSD({Q}, P, X) 3 
else 

X n ^:=X\{x n } 

Ho:=WRSD({Ti, ...,T n _i}, 5 do ,X„_ 1 ) 1 
G :=WRSD({Ti, . . . , T n _i}, 5 do ,X„_ 1 ) 2 
G:=G U map(t -> i U {T„}, G ) 
F:=F • WRSD({Ti, . . . ,T„_ 1 }, S do) X n _ 1 ) 3 
i:=0, S dl , +1 :=T n 
while ^ ^ do 

i:=« + l, H i: =0, G ; :=0 

Let Rdi be the d^-th principal subresultant coefficient of T n and P w.r.t. x n 
for H e H t i do 

H^I^UWRSDtH,^,^.!)! 
G. t :=G. t U WRSD(H, P rfi , X n _!) 2 
P:=F-WRSD(H,ii di ,X„_ 1 )3 

for G e Gi do 

H:=H U {G U {Sdi}}} 

Let Q be the pseudoquotient of T n and w.r.t. x n 
if deg(Q,x„) > then 

G:=G U WRSD(G U {Q}, P, X) 2 
P:=P • WRSD(G U {Q},P, X) 3 



41 return 



i,F 
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W^ttF, G}\I f ) = U^V^aS^R^. . . , R do }\I F R di ); 

(2)V im (F\Gl F ) = V-^jIFXGIfR^) U UV =1 V^ M ({Q (ii ,E (ii _ 1 , . . . , R do }\Gl F R di ). 

Proof. Assume that S^+i, S^, . . . , Si, So is the subresultant chain of F and G w.r.t. x n . Remark 
that S da = So = rcs(F, G), S^ — G, S^+i = F and S dv = Iq c G where c is a non-negative integer. 

(1) Without loss of generality, assume that V-^jjjj({F, G}\l F ) ^ 0. For any (ai, . . . ,a n ) £ 
Yjquj({F,G}\l F ), let b = (01, . . . , a n -i). If G(b) = 0, by the definition of principal subresul- 
tant coefficient, R di (b) = and thus R di {a) — for any i (1 < i < v). Hence, (oi, . . . ,a n ) £ 
^K\u]({ s d v+1 ,R dv ,- ■ -,Rd }\lFRd v+1 )- Ifdeg(G(6),x n ) > 0, since deg(G(6), x n ) < dcg(F(b),x n ) 
deg(F,x n ), it is reasonable to assume that the subresultant chain of F(b) and G(b) w.r.t. x n 
is S^+i, Sfx, . . . , S%, So and the associated principal coefficients are R^, . . . , Pi, Ro- Note 
that S^ = G(6),5 M +i = F(b) and by Lemma[21 we know that Sj(b) = l F (b) r: >Sj where rj is a 
non-negative integer for any j (I < j < /i + 1). According to the theories of subresultant chains, 
there exists an integer j (1 < j < fj,) such that Rj ^ and Rq = ■ . . = Rj-x = 0. Then Rj(b) ^ 
and Ro(b) = . . . = Rj-\(b) = 0. In addition, Sj is the greatest common divisor of F(b) and G(b) 
in X[C7J[a;„] and deg(Sj, x n ) = j. Hence Sj(a n ) = by F(b)(a n ) = G(b)(a n ) = 0. Note that 
deg(Sj,x„) = deg(Sj(b),x n ) = deg(Sj,x n ) = j, so there exists some i (1 < i < v) such that 
di =j. Therefore, (oi,...,a„) £Vj^({S di ,R di _ 1 ,...,R do }\I F R di ). 

On the other hand, for any (ai,...,o n ) 6 V-^jj-({5 d „ +1 , P<j„, • • • , Pd }\lFPd„ +1 ), let b = 
(ai,... , a n _i). As R di (b) = for any i (1 < i < u), G(6) = follows from Lemma [2 Hence, 
(ai,...,a n ) £ V^jjj{{F,G}\I f ). For any i (1 < i < v) and for any (ai,...,a„) £ V^jj-({5 di , 
R di _ 1 , ■ ■ ■ , R do }\I F R di ), it is not difficult to check (oi, . . . ,a„) G ({F, G}\Ip ) similarly as 

what has been discussed in the last paragraph. 

(2) The proof is similar to that of (I). □ 

Remark 4. If F and G are polynomials in K[X] (0 < deg(G,a;„) < deg(F,x n )), suppose 
S dv , . . . , S dl , S do is the regular subresultant chain of F and G w.r.t. x n . Let S dv+1 = F. Assume 
that R di is the di-th principal subresultant coefficient of F and G w.r.t. x n for any i (0 < i < v+l) 
and Q di is the pseudoquotient of F and S di w.r.t. x n for any i (1 < i < v). Similarly, we have 

(1) V({F, G}\l F ) = U^+ 1 1 V({S' di) i? rfi _ 1 , . . . ,R do }\I F R di ); 

(2) V(F\GI F ) = V(F\Gl F R do ) U uy =1 V({Q di , R di _ lt . . . ,R do }\Gl F R di ). 

Lemma 4. Let P £ K[U][X] and T = {Ti, . . . ,T n } be a zero-dimensional regular chain in 
K[U][X]. If S = rcs(PT) ^ 0, then res(P(o),T(a)) ^ for any a £ K d \V u (S tcs(1 t , T)). 

Proof. For any a £ K '\V^(5 res(I T , T)), we prove the conclusion by induction on n. When 
n = 1, it is not difficult to check that S'o(a) = ±Itj (a) fe res(P(a), Ti(a)) where k is a non- 
negative integer by the definition of resultant. Then res(P(a), T(a)) = res(P(a), 7i(a)) 7^ 
since So (a) 7^ 0. Assume that the conclusion holds when n < N (N > 1). If n = N, let 
Rn = res(P, T n ) and T„_i = {T u . . . , T^}. Since S = res(i?„, T„_i), res(Rn(a), T n _ x (o)) ^ 
by the induction hypothesis. Note that R n {a) = ±It„ (a) fc res(P(a), T n {a)) where k is a non- 
negative integer by the definition of resultant and T n _i(a) is a zero-dimensional regular chain in 
K[X] by LcmmaQ] By Remark[B;3), res(P(a), T(a)) = res(res(P(a),T„(a)),T n _i(a)) 7^ 0. □ 

Lemma 5. Given a zero-dimensional regular chain T = {Ti, . . . ,T„} m if[J7][X] anrf a poly- 
nomial P £ K[U][X], suppose P : = prcm(P,T). Then V^jr(T U {P}) = V^jj(T U {P x }) 
and V^t(T\P) = V^ J r(T\P 1 ). Furthermore, V(T(a) U {P(a)» = V(T(a) U {Pi (a)}) and 

V(T(o)\P(o)) = V(T(o)\Pi(a)) /or any a e ^\V a (res(I T , T)). 

Proof. According to the definition of successive pseudodivision, Pn™ =1 Iy. = ^™=i TiQi + Pi 
where Q±, . . . ,Q n Rre polynomials in I\ [U] a-nd ri , . . . , v n are non-negative integers. Note 

that for any a £ K d \ V u (res(Ix, T)), T specializes well at a by LemmaQ] Then it is not difficult 
to verify that Lemma [5] holds. □ 

Theorem 1. Algorithm^ terminates correctly. 
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Proof. The termination is similar as the termination of Algorithm RSD in [57J [35J . For a given 
zero-dimensional regular chain T = {Ti, . . . ,T„} in i4f[[/][X] and a polynomial P in i-C[[/][A], 
let WRSD(T, P, X) = H, G, F. Now we prove the correctness by induction on the recursive depth 
h of WRSD(T,P,X). 

When h = 1, according to Algorithm[TJ P is reduced w.r.t T and WRSD(T, P, X) can return at 
Line 5, Line 7, or Line 14. If WRSD(T, P, X) returns at Line 5 or Line 7, the conclusion follows from 
Lemma [TJ If WRSD(T, P, X) returns at Line 14, the conclusion follows from Lemma 01 Assume 
that the conclusion holds when h < N (N > 1). Suppose h = N. Then WRSD(T, P, X) can 
return at Line 3, Line 12, or Line 41. If WRSD(T,P, X) returns at Line 3, the conclusion follows 
from the induction hypothesis and Lemma[5j If WRSD(T, P, X) returns at Line 12, the conclusion 
follows from the induction hypothesis. Now we prove the conclusion when WRSD(T, P, X) returns 
at Line 41, which means P is reduced w.r.t. T, mvar(P) = x n and res(P, T) = 0. Suppose 
S^+i, Sp, . . . , Si, So is the subresultant chain of T n and P w.r.t. x n in K[U, X n -i] [x n ] where 
X n -\ = X\{x n } and Sd v , . . . , S^, Sd is the associated regular subresultant chain. Note that 
deg(P„,a;„) > deg(P,x„) > since P is reduced w.r.t. T. 

If n = 1, Srfj is the greatest common divisor of Ti and P in ii"[i7][xi] and hence V^r^ ({Ti , P}) 
= V^ryj({iS'd 1 }). Then the claim (1) in the specification of Algorithm Q] holds. Suppose Q is 
the pscudoquoticnt of T n and S dl w.r.t. xi. Remark that deg(Q, xi) > and there exists a 
positive integer k such that k > 2 and l| Ti = Sd^Q- Thus V^pr(T) = Vj^jjj(Sd 1 Q). Note 
that V^(5 dl ) C V^P). So V^(T\P) = V^S^QV) = V^(Q\P). Therefore 
the claim (2) in the specification of Algorithm [T] follows from the induction hypothesis. Re- 
mark that ls dl is a factor of Iq and according to Algorithm [TJ Iq is a factor of P. Thus for 

any a G K d \V u (F), Iti (a) ^ and deg(P(a),xi) > di > by the definition of subresultant. 
According to Lemma [21 S dl (a) is the great common divisor of P(a) and Ti(a) in iffzi]. Thus 
V(T(o) UP(o)) = V(S dl (a)). Since l Sdl (a) k Ti(a) = S dl (a)Q(a) and V(S dl {a)) C V(P(a)), 
V(T(a)\P(a)) = V{S dl (a)Q(a)\P(a)) = V(Q(a)\P(a)). Therefore the claim (3) in the specifi- 
cation of Algorithm [1] follows from the induction hypothesis. 

If n > 1, let Sd v+1 = S^+i and T„_i = {Ti, . . . , T„_i}. Suppose P^ is the principal 
subresultant coefficient of T n and P w.r.t. x n for any i (0 < i < v + 1) and assume that 
H = WRSD(T n _i,Sd , X n _i)i and G = WRSD(T„_i, S do , X n -i) 2 . Remark that H ^ because 
rcs(P do ,T n _i) = rcs(P,T) = 0. For any i (1 < i), let H* = U HeHl _ 1 WRSD(H, R dz , X„_i)i and 
G, = UHeHi_iWRSD(H, P^ , X n ^i)2 until there exists an integer / (1 < Z < u+l) such that H; = 0. 
That means H; = and Mj ^ for any j (0 < j < I). We can always get this integer I owing 
to the fact that Sd v+1 = T n . Then we have two sequences Ho, Hi, . . . ,Hj and Go, Gi, . . . , G;. 
Let Li = {i\l < i < 2,G, ^ 0}. According to AlgorithmEQ the first output of WRSD(T, P, X) is 
H = UigLj UGgG; (G U {Sdi})- It is not difficult to see that H is a finite set of zero-dimensional 
regular chains in A'[P][A]. By Lemma GUI), we know that 

V W (TU{P}) 

= V 7?M (T„_ 1 ) n V^(T„ U {P}\IrJ 

= V 7qUT( T «- 1 ) U i=l V KpT^' S '*'" R *- 1 '' ' -'^olVT^PdJ 

= u^ 1 1 (v 7 ^(t„_i) n v 7 ^({s , di ,p (ii _ 1 , . . . ,R do }\l Tn Rdi)) 

For any i (1 < i < u + 1), If i g Li, according to the induction hypothesis and the construction 
of Gi, we get 

V ifp]"( T «-i) n ^ii[U]({ s d i ,Rd^ 1 ,- ■ ■ , R do }\lT„Rdi) 
= V kT77]"( T "- 1 ) n ^^Ujd^ditRdi-!, ■ ■ ■ >Rdo}\Rdi) 
= v 7fTT7r( Tri - 1 u i Rd «}) n ^K\u]{{ s d i ,Rd^ 1 ,- ■ .,R dl }\R di ) 
= U He H V^pj-(H) n V^u]{{Sdi, Pdi_!, • ■ . ,Pd 1 }\Pd i ) 
= UHeH V^pj-(H U {Pdj}) n V^yj-({S'd i ,Pd i _ 1 , • • • , Rd 2 }\Rdi) 

= UHeH!V^jr(H) ("1 y^Tfw({S di , Rdi-i, ■ ■ ■ ,Rd 2 }\Rdi) 
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= . • . = Uhgh^V^H) n Vj^iiSd, }\R dl ) 

= u HeHl _ 1 v 1?M (H\i? di ) n Vj^dSd,}) 

= U G6 g 4 V^(G) n V-^dSd,}) 
= UGec i V IM (Gu{S dj }). 

If I < i < v + 1. according to the induction hypothesis and = 0, similarly, wc know that 

V KpI^ T "~ 1 - ) n ^i^iujiiSdiiRdi-!, ■ • • , Rdo}\lT n Rdi) 
= Uhgh,V 7 ^(H) n VjcffliiSditRdi-!,- ■ ■ , Rdi +1 }\Rdi) 
= 0. 

If 1 < i < I and i £ L\, similarly, we get 

= u Hg H 4 _ 1 V^ T (H) n VjcmdS^Rd,) 
= UHeH^xV^HViidj n V^iSdJ) 
= Uoe^v^G) nv^s*}) 

= 0. 

Therefore, V 7? pj(T U {P}) = U ieLl U GeGl V^(G U {S di }) = U HeH V 5? p ] -(H) and hence the 
claim (1) in the specification of Algorithm [1] holds. Furthermore, as discussed above, we figure 
out that V 1? pj(TU{P}) = if and only if H = 0. Actually, when res(P,T) = 0, V-^jj (TU {P}) 
can't be according to Proposition [^3) and thus H ^ 0. Similarly, we can prove that the claim 
(2) in the specification of Algorithm [1] holds on the basis of Lemma [3J2). Besides, it also can be 
shown that V^T\P) = if and only if G = 0. 

For any aff \V U {F), it is easy to check that H specializes well at a for any H £ I by 
the induction hypothesis and we only need to prove that V(T(a) U {P(a)}) = UHenV(H(a)). 
If deg(P(a), x n ) = 0, by the definition of principal subrcsultant coefficients, Rd i (a) = for 
any i (1 < i < v). Thus it is not difficult to check that Gj = 0(l<i<i;) and hence 
H = U GeG „ +1 {G U Note that R do (a) = I Tn ( )deg(i J > *»)p( )deg(:r„,x n ) ) so V (T(a) U 

{P(a}}) = V(T(a) U {R do (a)}). By G* = (1 < i < v) and the induction hypothesis, 

U GeGll+1 V(G(a)U{S dv+1 (a)}) 

= U He iJ(H( tt )\JJ d „ +1 («))nV(T„(fl)) 

= U H eH„_ 1 V(H(a)) n V({r n (a), J R d „(a)}\ii d „ +1 (a)) 

= U He H„_xV(H(a)) nV(T„(a)\i? dt;+1 (a)) 

= ■ • ■ = U HG H V(H(a)) n V({T„(o), i? dl (a)}\i? dt , +1 (a)) 

= U He H V(H(a)) nV(T n (a)\R dv+1 (a)) 

= V(T n _i (a))) n V({T„(a), i? do (a)}\i? d „ +1 (a)) 

= V(T(a)U{ J R do (a)} 

Therefore, V(T(a) U {P(a)» = U HeH V(H(fl)). 

If deg(P(a), cc„) > 0, it is reasonable to assume that the subresultant chain of T n {a) and 
P(a) w.r.t. x n is S^+i, S^, . . . , So- By Lcmma[2J we know that Si(a) = I<r n (a) 7 * Si where Ti is 
a non-negative integer for any i (0 < i < fi + 1). Suppose L2 = {i|l < i < u + 1, R di (a) 7^ 0}. 
It is not difficulty to check that L\ C L2 by the induction hypothesis and it is reasonable to 
assume that L2 = {ji, ■ ■ ■ ,jk,jk+i} (k > 1) such that < dj 1 < . . . < dj k < dj k+1 = d v+ i. Then 
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So, Sd Jl , ■ • ■ , Sd jk is the regular subresultant chain of T n (a) and P(a) w.r.t. x n . By Rcmark^l), 

V(T(a)U{P(a)}) 
= V(T n _i(o)) D V(T n (o) U {P(o)» 
= V(T B _i(o)) n V(T n (o) U {P(o)}\I Tb(o) ) 
= V(T B _i(o)) n U^VtfS^ , Rd^ R do }\l Tn(a) R dh ) 
= U^(V(T n ^(a)) n V({5 djt , i? djt i , . . . , R do }\l Tn{a) R djt )) 
= U it6i3 (V(T n _x(a)) n V({Sd H , , . . . , R do }\l Tn(a) Rd H )) 
For any j t £ L2, if Jt £ Ti, then by the induction hypothesis, 

V(T n _l(o)) fl V({S djt , , . . . , P^ , i? do }\lT„(a)Pd 3t ) 

= V(T„_i(a)) n V({5 djt (a), R d]t _ 1 (a), R d]1 (a), R do (a)}\R d]t (a)) 
= UHgH^VCHO*)) n V({S Jt (a)}\RdJa)) 
= U He M, t _ 1 V(H(a)\P dit (a)) n V({5 djt (a)}) 
= U GeGj( V(G(a))nV({S dji (a)}) 
= U GeG3( V(G(a)U{5 djt (a)}) 

For any j t £ L^\Li, if jt < Z, by the induction hypothesis, UHgH _ 1 V(H.(a)\R d - (a)) = since 
% = 0. Then 

V(T n _i(o)) n V({S du , R d]t _ t , • . • , P do }\lT„(a)-Rd, t ) 

= V(T„_ x (a)) n V({5 djt (a), P djti (a), . . . , P do (a)}\P djt (a)) 
= UHGHfc.xVWo)) n V({5 d , t (a)}\i? djt (a)) 
= Uh^,., V(H(a)\i? d . t (a)) n V({5 djt (a)}) 
= 0. 

If jt > I, by the induction hypothesis, UhgH;_i V(H(a) U {R dl (a)}) = since = 0. Then 

V(T n _!(a)) n V({Sd H , Rd^ R do }\I Tn(a) R dji ) 
= V(T n _!(o)) n V({S dH (a), P^ (a), . . . , R do (a)}\R dh (a)) 
= Uhei,., V(H(a)) n V({5 d . ( (a), P^ (a), . . . , P d;+1 (a), R dl (a)}\R djt (a)) 
= 

Therefore, V(T(o) U {P(a)» = U A6£l U GeGjf V(G(a) U {S»(o)» = U He HV(H(a)) and hence 
the claim(3a) in the specification of Algorithm [1] holds. It is similar to check that the claim (3b) 
in the specification of Algorithm [T] holds by Remark|21j2). □ 

3.2 Converting To Regular Chains 

Let T be a triangular set in P[[/][A]. If mvar(T)=A, we can compute a finite set G of zero- 
dimensional regular chains in P[Z7][A] and a polynomial F £ K[U] on the basis of Algorithm 

□ such that V 7?m (T\It) = U GeG V^ 7r (G) and for any a £ K d \V u (F), V(T(o)\I T (o)) = 
UgggV(G(cs)) and G specializes well at a for any G £ G if G ^ 0. The algorithm is presented 
as Algorithm [21 which plays a key role in Algorithm [3] proposed in the next section. 

Theorem 2. Algorithm^ terminates correctly. 

Proof. If the input T is a regular chain, then the termination holds obviously and the correctness 
follows from Lemma [TJ Now we assume that T is not a regular chain and let k be the minimal 
integer k (1 < k < n) such that T& = {T±, . . . , T^} is a regular chain and {Ti, . . . , Tk+i} is not a 
regular chain. Remark that this assumption is reasonable owing to the fact that at least {Ti } is 
a regular chain in K[U}[X\. Let T„_ fe = {T k+1 , . . . ,T„}. 
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First of all, we prove the termination. Assume that ZDToRC(T, A) doesn't terminate. Then 
we can get at least one regular chain R e WRSD(Tfc, lr k+1 , {xi, ■ . ■ , x k })2 such that ZDToRC(R U 
T„_k,X) can't terminate. According to Algorithm^ there exists k.2 (1 < k^ < n) such that fc 2 
is the minimal integer such that RU {T k +i, . . . , T k2 } is a regular chain but RU {T k +i, ■ ■ . , Tk 2 +\] 
is not a regular chain. It is easy to check that fc 2 > k. Since ZDToRC(R U T n -k,X) doesn't 
terminate, the rest can be done in the same manner and we can get an infinite sequence of 
positive integers k = ki < fc 2 < ■ • ■ < k t < .... Note that all positive integers in this infinite 
sequence must equal to or be less than n and it's impossible obviously. Therefore, Algorithm^ 
terminates. 

Now we prove the correctness by induction on the recursive depth h. If h = 1, which im- 
plies WRSD(T fc ,I Tfc+1 ,{a;i,...,x fc }) 2 = 0, then ZDToRC(T,A) = 0, WRSD(T fc , l Tk+1 , {x u . . . , x k }) 3 . 
Let F = ZDToRC(T, A) 2 . According to Algorithm [TJ V^pjr(T k \I Tfe+1 ) = and for any a 6 

K d \V u {F), V(T k (a)\I Tfe+1 (a)) = . That means Vj^{T\l T ) = and for any a e R d \V u (F), 
V(T(a)\I T (o)) = 0. Thus the specification of Algorithm [2] holds. 

Assume that the specification of Algorithm [5] holds for h < N (N > 1). If ft. = N, assume 
that ZDToRC(T,A) = G,F. Let W=WRSD(T fe , Ir fc+I , {xi, ■ . ■ ,x k }). Then we can assume that 
W2 is a non-empty set {Hi,...,H„i} of zero-dimensional regular chains in K[U][xi, . . . , Xk\- 
Let Ri = Hi U T„_fc for any i (1 < i < rri). Note that W3 is a factor of F. By Algo- 

rithmm V^TAItvJ = U™^^^) and for any a € K d \V u (F), V(T k (a)\l Tk+1 (a)) = 
U^VCH^a)). Let L = {i\l < i < m, V^^RAIrJ + 0}- Then G = U ieL ZDToRC(Ri, X) x and 
ZDToRC(Ri, X)<z is a factor of F for any i (1 < i < m). Thus by the induction hypothesis, we 
have 

%[o](T\It) 

= v 1?M (T fc _ 1 \i Tfe+1 ) n v w (T n _ fc \i Tn _J 
= u™ 1 v im {n i ) n v m (T„_ fe \i T „_J 
= u™ 1 v 1?m (rair i ) 

= U^lV^RAIrJ 
= U G eGV^jj]-(G). 

Hence, the claim (1) in the specification of Algorithm [5] holds and from the proof one can see 
that V^r(T\I T ) = if and only if G = 0. If G ^ 0, every G in G must be returned by Line 2 
in some recursion according to Algorithm [2] For this reason, res(Ic G) is a factor of F. Hence 
for any a £ K d \V u (F), G specializes well at a by Lemma [TJ Furthermore, whether G = or 
not, by the induction hypothesis we get 

V(T(o)\I T (o)) 

= V(T fc _i(o)\I^ +1 (o))nV(T n _ fc (o)\I TB _ k (o)) 
= U™ 1 V(H l (a))nV(T„_ fe (a)\lT„_ fc (a)) 
= U™ 1 V(R l (a)\I Rl (a)) 
= U < eiV(R i (o)\lR 4 (a)) 
= U GeG V(G(a)). 

Hence, the claim (2) in the specification of Algorithm [5] holds. □ 

3.3 Computing RDU Varieties 

We present the main result in this section. Algorithm [3] shows how to compute a generic 
regular decomposition and the associated RDU variety of a given generic zero-dimensional system 
simultaneously. As a matter of fact, whether a given system is generic zero-dimensional can be 
checked by the Algorithm [3] itself. For this reason, we just assume that the input system of 
Algorithm [3] is always generic zero-dimensional. 
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Algorithm 2. ZDToRC 



Input: A triangular set T = {Ti, . . . , T„} in A[L/][A] satisfying mvar(T) = X, variables 
X = {x%, . . .,x n }. 

Output: A finite set G of zero-dimensional regular chains in A[t/][A] and a polynomial 
F e K[U], such that 

(1) V^(T\I T ) = Ug^V-^G); 

(2) for any a G K d \V u (F), V(T(a)\I T (a)) = U GeG V(G(a)) and G specializes well at a for 
any G G G if G ^ 0. 

1 if T is a regular chain then 

2 | return {T}, res (I T ,T) 

3 else 

Find the minimal integer k (1 < k < n) such that T& = {Ti, . . . , Tfc} is a regular chain 
and Tk+i = {Ti, . . . , T&, Tt + i} is not a regular chain 

5 W: =WRSD(T /£ ,I Tfc+1 ,{x 1 ,...,x fc }) 

6 if W2 = then 

7 I return 0, W3 

8 F: =VF 3 , G: =0 

9 for T in W2 do 



10 R: ={op(T),T fe+1 ,...,T n } 

11 |_ G: =GUZDToRC(R, A)i, F: =F-ZDToRC(R, X) 2 

12 return G, F 



Lemma 6. Given a zero- dimensional regular chain T in K\U][X\ and a finite set of poly- 
nomials P C K[U][X], suppose V^pjj-(T) C V-^(P). Then V(T(a)) C V(P(a)) /or any 

ae jT\v^(res(I T) T)). 

Troo/. By Proposition E(3), V W (T\I T ) - V^(T). Then V 7?M (T\I T ) C V^P). By 
Proposition [2^2) , P C sqrt(sat(T) K j ;7 ]j X j) and hence for any P G P, there exist a positive 
integer A: such that P k G sat(T) K j (7 jj X j. By Proposition H]T), prem(P fe ,T) = 0. Remark that 
P fc G A'[C/][A] and T C K[U][X], so V^(T\I T ) C Vj^P*) = V^(P). For any a = (a u . . . , a d ) G 
A \V c/ (i'cs(It, T)), T(a) is a zero-dimensional regular chain in K[X] and It(o) 7^ by Lemma 
[TJ Hence res(I T (o), T(o)) = rcs(I T(a) , T(a)) ^ and thus for any b = (6 X , . . . , & n ) G V(T(a)), 6 £ 
V(I T (a)). That implies (ai, . . . , a d , h, . . . , b n ) G V^(T\I T ) C V^-(P). Therefore, 6 G V(P(a)) 
and then V(T(a)) C V(P(a)). □ 

Theorem 3. Algorithmic terminates correctly. 

Proof. Since the termination follows from the termination of Algorithm [T] and Algorithm [2j 
we only need to prove the correctness. Assume that P is a generic zero-dimensional system in 
K[U] [X] and {Ci, . . . , C m } is a Wu's decomposition of P C A[[/][A] computed by Wu's method. 
According to Wu's method and Algorithm^ we know that the claim (1) in the specification of 
Algorithm [3] holds. 

Now we prove the claim (2) in the specification of Algorithm [3] by induction on m. If m = 1, 
Ci is a characteristic set of P in A[C/][A] and Ci = {C\} C K[U] by Wu's method. That means 
V-^jj-(P) = 0. According to Algorithm |3J the first output of RDUForZD(P, X) is and the second 

output is exactly C\. In fact, for any a £ V u (Ci), C\{a) G K. Note that C\{a) G (P(o))k[X1 
by Ci G (P)k[u][x}- Thus V(P(a)) C V(Ci(a)) = and the claim (2) in the specification of 
Algorithm |3j holds. 

Assume that the conclusion holds for m < N (N > 1). If m = N, suppose Ci is the 
characteristic set of P computed by Wu's method. Since m > 1, we know that Ci <£_ K[U] 
and VjfptfP) = V^C^IcO U V^P U d U {I Cl }). Let RDUForZD(P, X) = T, F and 
ZDToRC(Ci, X) = Ti,Fi. Then Tj C T and F-y is a factor of F. In addition, if RDUForZD(P U 
Ci U {I C J, X) = T 2 , F 2 , we have T 1 \JT 2 = T and F = F x ■ F 2 by Wu's method. We only prove 
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the conclusion when Ti 7^ and T2 7^ 0. In fact, if Ti =0 or T2 = 0, the proof is similar. 
For any a £ V U (F), a £ V U (F 1 ) and then V((P U Ci U {I Cl })(a)) = U TG T 2 V(T(a)) by the 
induction hypothesis. Therefore, in order to prove V(P(a)) = UxeTV(T(a)), we only need to 
show V(P(a)) = U T£ T 1 V(T(a)) U V((P U Ci U {I Cl })( a ))- As a matter of fact, by Algorithm 
m U TGTl V(T(a)) = V(C 1 (a)\I C i(a)). Note that d(a) G (P(a))^ [X] , so V(P(o)) C V(C 1 (a)). 
Then V(P(a)) C V(Ci(a)\Ioi(o)) U V((P U^U {I Cl })(a)) = U TeTl V(T(a)) U V((P UCjU 
{Ici})(a))> On the other hand, by the claim (1), V ^.^ (P) = UtetV ^^j (T) and thus for any 
T G Ti C T, V^pr(T) C V-^r(P). According to Algorithmic we know that res(I T , T)(a) ^ 0. 
By LcmmaEl Uts^ V(T(a)) U V((P U Ci U {I Cl })(a)) C V(P(a)) and we arc done. □ 



Algorithm 3. RDUForZD 



Input: A generic zero-dimensional system P in A'[[/][A], variables X = {xx, . . . , x n }. 
Output: A set T of zero-dimensional regular chains in if[£/][A] and a polynomial 
F C K[U], such that 

(1) V M (P) = UtetV^CT); 

(2) for any a £ K d \V u (F), V(P(a)) = U Te TV(T(a)) and T specializes well at a for any 
T 6 T if T / 0. 

1 Compute a Wu's decomposition § = {Ci, . . . , C rn } of P in if[t/][A] by Wu's method 

2 § 2 :=§\{C m }, F:=op(C m ) 

3 if § 2 = then 

4 |_ return 0,F 

5 T:=0 

6 for C in S2 do 
7 

8 



W:=ZDToRC(C,A) 
T:=TUWi, F:=F-W 2 



9 return T, F 



For a generic zero-dimensional system P in if [C/][X]. Suppose Vj^jy(P) 7^ 0. Assume that 
§ = {Ci, . . . , C m } is a Wu's decomposition of P C -ftT[£/][A] computed by Wu's method and let 
W = {Ce MlVj^CVc) + 0}. Remark that W/J since Vj^P) ^ 0. Obviously, W is a 
non-redundant parametric triangular decomposition of P in if[[/][X]. If F = RDUForZD(P, A~) 2 , 
we have the following corollary. 

Corollary 1. The non-redundant parametric triangular decomposition W of P in K[U][X] is 
stable at a if a G K d \V u (F). 

Proof. Assume that W = {Si, . . . , SJ. For any i (1 < i < t), assume that T, = ZDToRC(S,, X)i 
and according to Algorithm^ we know that ZDToRC(Si, A)2 is a factor of F. Then 7^ 0, since 
V^7yr(Si\Is i ) 7^ 0- Firstly, we need to interpret that for any a V U (F), Is;(a) 7^ 0. In fact, 
if l Si (a) = 0, V(S l (a)\I Sl (a)) = 0. On the other hand, by Algorithm [5] and Algorithm El we 
know that V(S 4 (a)\I Si (a)) = U T eT 4 V(T(a)). Note that T.; ^ and T(a) is a regular chain for 
any T G Tj, so U T eT s V(T(a)) ^ by Proposition [TJ That implies V(S.;(a)\I Sl (a)) ^ 0, a con- 
tradiction. Therefore, by Algorithm^ V(P(a)) = U* =1 Utgt, V(T(a)) = Uf =1 V(S;(a)\I Si (a)) = 
U* =1 V(Si(a)\Is i ( a )) and hence the conclusion holds. □ 

Corollary [1] indicates that we can also obtain a non-redundant parametric triangular de- 
composition of a given generic zero-dimensional system and the decomposition is stable at any 
parameter value that is not on the RDU variety computed by Algorithm [3] The following Ex- 
ample Q] is presented to illustrate how Algorithm [3] and Corollary Q] work out. 

Example 1. Consider the system 

!(u - l)x\ + (x\ - 2uxi + u 2 + l)x 2 + x\ — Xl, 
(xi-u)(x 2 + l), 
{xx - uf 
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where x± and xi are variables (x± -< X2) and u is a parameter. 

Step 1: According to Wu's method, we compute a Wu's decomposition § = {Ci, C2, C3} of 
P in R[u][xi, X2] where Ci = {{x\ — u) 2 , [x\ — u)(x2 + 1)}, C 2 = {x\ — u, (u — \ )x\ + X2 + u 2 — u) 
and C 3 = {u — 1}. 

Step 2: Let § 2 = {Ci, C 2 } and F = op(C 3 ) =u-l. 

Step 3: Because § 2 ^ 0, we need to execute ZDToRC(Ci, {x 1} x 2 }) and ZDToRC(C 2 , {xi,x 2 })- 

Step 3.1: It is easy to see that Ci is not a regular chain in M.[u] [xi, x 2 ]. By calling 
WRSD({(a;i -u) 2 },xi -u, {xi}), we obtain H = {{xi -u}}, G = and F 1 = 1. Then the first and 
second outputs of ZDToRC(Ci, {x\ 1 x 2 }) are and 1 respectively. That means Vjj^(Ci\xi — u) = 
and for any a £ C, V(d(a)\xi - a) = 0. 

Step 3.2: The first and second outputs of ZDToRC(C 2 , {a;i,x 2 }) are {C 2 } and u — 1 
respectively owing to the fact that C 2 is a regular chain in a; 2 ]. In addition, we know 

that for any a £ C, C 2 (a) specializes well at a as long as a ^ 1 by Lemma [T] 

Combining the above results together, we get a generic regular decomposition {C 2 } of P 
in R[u][xi,a; 2 ] and a RDU variety V = {a £ C|o —1 = 0} such that for any a £ C\V, 
V(P(a)) = V(C 2 (a)) where C 2 (a) is a regular chain in C[xi,x 2 ]. That means for almost all 
parameter values, we can express the solutions of P by C 2 . Furthermore, It is should be noted 
that Algorithm [3] eliminates a redundant branch Ci from § 2 and this result should be attributed 
to weakly relatively simplicial decomposition. As a result, we also get a non-redundant para- 
metric triangular decomposition {C 2 } of P. As what Corollary [T] shows, this non-redundant 
decomposition is stable at any a £ C\V. 

4 Implementation 

We have implemented Algorithm [3] as a function rduforzd on the basis of DISCOVERER 
P5] using Maple 15. More specifically, Wu's method for computing parametric triangular de- 
compositions introduced in Section [2] is implemented as a function wusolve and Algorithm [1] is 
implemented as a function wrsd. Throughout this section, all the results are obtained in Maple 
15 using an Intcl(R) Core(TM) i5 processor(3. 20GHz CPU and 4GB total memory). 



Table 1. Comparing wrsd and rsd 



system 


U 


X 


time 


H 


G 


F 


wrsd 


rsd 


wrsd 


rsd 


wrsd 


rsd 


wrsd 


rsd 


EX1 


3 


7 


0.203 


0.296 


1 


1 


2 


2 


11 


11 


EX2 


2 


6 


0.156 


0.234 


1 


1 


2 


2 


9 


9 


EX3 


4 


6 


0.437 


1.061 


1 


1 


2 


2 


10 


10 


EX4 


3 


6 


1.482 


1.452 


1 


1 


1 


1 


11 


11 


EX5 


4 


4 


0.281 


0.281 


1 


1 


1 


1 


9 


9 


EX6 


4 


4 


0.218 


0.218 


1 


1 


1 


1 


9 


7 


EX7 


3 


3 


0.327 


0.375 


1 


1 


1 


1 


8 


8 


EX8 


3 


3 


0.468 


0.577 


1 


1 


1 


1 


7 


7 



We run 8 examples^ using wrsd and rsd on the same computer with Maple 15 and the 
comparisons about timings and results are presented in Tabic [U where columns X and U represent 
the cardinal numbers of the variables and the parameters, respectively, column time reports the 
timings in seconds, columns H and G represent the numbers of branches in the first and second 
outputs, respectively, column F represents the numbers of irreducible factors over the field of 
rational numbers of the third output. The empirical data shows that wrsd performs as well as 
rsd with higher efficiency in most cases. 

We also run several examples from the literature using rduforzd with Maple 15 and part of 
the empirical data about timings is presented in Tabled Therein, systems 1-4 are selected from 
[T3"] . systems 5-6 and 12-14 are provided by Sun [TU], systems 7-11 are modified from Examples 

4 http: / /www. is. pku.edu.cn/~xbc/ExForRSD.txt 

5 Algorithm RSD in [27] was implemented as a subfunction rsd in DISCOVERER by Xia [23] , 
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2-4 in [5] and systems 15-22 are provided by Chen [3]. In Tabled columns X and U represent 
the cardinal numbers of the variables and the parameters, respectively, column wusolve reports 
the timings in seconds cost by computing Wu's decompositions, column zdtorc represents the 
timings in seconds cost by computing weakly relatively simplicial decompositions and some other 
steps required in Algorithm 03 and column total shows the timings added by the timings in the 
former two columns. It is indicated that our method can be applied to a wide range of practical 
problems with reasonable time cost. Furthermore, it is interesting to note that computing generic 
regular decompositions and the associated RDU varieties do not require much more time cost 
than Wu's decompositions when solving practical problems as shown in Tabled 



Table 2. timings o/rduforzd 



number 


system 


U 


X 


wusolve 


zdtorc 


total 


1. 


S5 


4 


4 


0.016 


0. 


0.016 


2. 


S9 


3 


3 


0.047 


0. 


0.047 


3. 


S16 


3 


12 


0.031 


0. 


0.031 


1 


S18 


3 


2 


0.437 


0.046 


0.483 


5. 


SY1 


3 


2 


0.016 


0.015 


0.031 


6. 


SY2 


4 


1 


0.015 


0. 


0.015 


7. 


SCC1 


3 


4 


0.016 


0. 


0.016 


8. 


SCC2 


4 


7 


0.015 


0.016 


0.031 


9. 


SCC3 


6 


11 


0.110 


0.078 


0.188 


10. 


SCC4 


4 


7 


0.046 


0.032 


0.078 


11. 


SCC5 


4 


5 


0.032 


0.015 


0.047 


12. 


P3P 


5 


2 


0.031 


0.016 


0.047 


13. 


F4 


4 


1 


0.016 


0. 


0.016 


11 


F6 


4 


1 


0.015 


0. 


0.015 


15. 


Gcrdt 


3 


4 


0.031 


0. 


0.031 


16. 


Wang93 


2 


3 


0.062 


0. 


0.062 


17. 


Ley kin- 1 


4 


4 


0.078 


0. 


0.078 


18. 


Neural 


l 


3 


0.078 


0.016 


0.094 


19. 


Pavelle 


4 


4 


0.140 


0.094 


0.234 


20. 


gcnLinSyst-3-3 


12 


3 


0.032 


0. 


0.032 


21. 


AlkashiSinus 


3 


6 


0.015 


0. 


0.015 


22. 


LanconeLLi 


7 


4 


0.093 


0. 


0.093 



5 Conclusions 

We give an algorithm for computing generic regular decompositions and the associated RDU 
varieties simultaneously for generic zero-dimensional systems in this paper. That means for 
almost all parameter values, we can judge whether a given generic zero-dimensional system 
generally has solutions and present the solutions in an algebraic closed field by finitely many 
zero-dimensional regular chains if it does. As a result, questions (1) and (2) in Section Q] arc 
answered to some extent. 

Remark that although we refer to many methods for computing triangular decompositions 
such as Wu's decomposition, relatively simplicial decomposition and regular chain decomposition 
in our discussion, our main direction is not to give a new method for computing triangular 
decompositions but to provide an idea to discuss the relationship between parameter values 
and representations of solutions for parametric systems. Following the idea presented in this 
paper, any algorithm for computing regular chain decompositions can be probably adopted to 
computing generic regular decompositions and the associated RDU varieties if the algorithm is 
based on computing resultants and pseudorcmaindcrs. Algorithm [3] can be modified so that it 
can be applied to deal with general systems besides generic zero-dimensional ones. Furthermore, 
we can answer questions (1) and (2) step by step in algebraic closed fields completely by using 
Algorithm^ Actually, the concept of RDU variety is closely bound up with the concept of border 
polynomial (BP) in [24l [25j [26]. More precisely, if V is a RDU variety of P w.r.t. some generic 
regular decomposition, there exists a related border polynomial F of P such that V is contained 
in the hypersurface generated by F. A clearer characterization of the relationship between BPs 
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and RDU varieties is interesting but not what we concern here. We'll discuss these topics in the 
future. 
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